Cavity nonlinear optics with few photons and ultracold quantum particles 
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The light force on particles trapped in the field of a high-Q cavity mode depends on the quan- 
tum state of field and particle. Different photon numbers generate different optical potentials and 
different motional states induce different field evolution. Even for weak saturation and linear polar- 
izability the induced particle motion leads to nonlinear field dynamics. We derive a corresponding 
effective field Hamiltonian containing all the powers of the photon number operator, which predicts 
nonlinear phase shifts and squeezing even at the few-photon level. Wave-function simulations of the 
full particle-field dynamics confirm this and show significant particle-field entanglement in addition. 



It is commonly assumed that ofFresonant interaction 
of coherent light with a linear polarizable medium only 
creates linear phase shifts and coherent states of the scat- 
tered light. This is especially the case for weak light fields 
involving only few photons, where generating any nonlin- 
ear phase shift is a serious challenge [l| while such devices 
would be highly desirable in photonics or quantum in- 
formation applications 0]- One possibility to generate 
nonlinear phase shifts even at the single photon level was 
demonstrated via resonant strong coupling of an atomic 
transition and a photon in a high- Q cavity Q . Recently 
it was noted that nonclassical light fields can also emerge 
from linear scattering off a weakly excited medium, if its 
motional state has genuine quantum properties [4;] . Here 
we combine those two ideas and study the effective non- 
linear cavity field dynamics generated from the photon 
number dependence of light forces on the quantum par- 
ticle motion in the cavity-enhanced optical potential [5|. 

In a prototype setup an ultracold particle is placed in 
the optical potential generated by a cavity field mode. 
For large detuning between atomic and cavity resonance 
spontaneous emission is small and the cavity field simply 
generates a photon-number dependent trapping poten- 
tial. For any superposition state involving several pho- 
ton numbers, as e.g. a coherent state, the localization of 
the particle thus differs for each photon number. As the 
particle position spread in turn determines effective re- 
fractive index and field phase shift this implies effective 
nonlinear field evolution. Experimentally a closely re- 
lated nonlinear effect was found in a pioneering setup in 
Berkeley |6|] . As several similar setups couplinga EEC to 
a high-Q cavity mode have just been realized [7[, further 
detailed experimental studies will emerge soon. 

Naturally the coupled particle-field dynamics gener- 
ates entanglement, which partly persists even in a steady 
state, where each photon number is correlated with a dif- 
ferent particle state. Interestingly, even after tracing over 
the particle degrees of freedom the field exhibits nonclas- 
sical properties and its dynamics can be approximated by 
a simple effective field Hamiltonian. We study its central 
predictions, and check its validity by a full particle-field 



FIG. f : (Color online) Scheme of the system with a quan- 
tum particle in the field of a driven cavity. A single driven 
standing wave mode generates optical potentials of different 
depth depending on the photon number. We also display the 
corresponding particle ground-state wavefunctions. 



wave-function simulation. 

Our prototype CQED system is depicted in Fig. [1] 
Following standard assumptions for a particle moving 
in ID along the axis of a far red detuned single cavity 
mode, we start from a Jaynes-Cummings Hamiltonian 
with eliminated internal state(s) [1, Q to get {fi= 1): 
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Here x, p, and are the particle position and momen- 
tum operators and its mass; a is the field operator, K 
is the mode wave number, Ac = uj ~ ujc is the cavity 
detuning (w is the laser, wc is the mode frequency), rj 
is cavity pump amplitude. Vb describes an additional 
c-number potential along the cavity axis, it ensures a 
particle bound state in the absence of cavity photons, 
but can be neglected for larger photon numbers. A cen- 
tral parameter here is Uq, the optical potential depth per 
photon, which also gives the cavity frequency shift per 
article and is proportional to the particle susceptibility 
As we consider small cavities, photon loss is impor- 
tant even for good mirrors and will be described by the 

standard Liouvillean [lo| Lp — k (^2apa^ — [a^a, p]^^, 

where k is the mode linewidth. In the following dJl (H) 
stands for a Master equation defined by Hamiltonian H 
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and this Liouvillean. 

The most interesting part of the dynamics arises from 
the second term of the Hamiltonian ([T]) describing dis- 
persive particle-field interaction. Its effect can clearly be 
seen from an expansion of the wavefunction in a product 
basis l^*) — J2n'^n\n) \(t)n), where \n) is the n-photon 
field state and a particle wavefunction, which will 
evolve subject to the potential V(n) f{Kx) \\^. Natu- 
rally, this creates entanglement between the particle and 
field on the timescale of the particle motion, character- 
ized by Wrcc- The pump term proportional to r] in the 
Hamiltonian ([1]) and the photon loss described by quan- 
tum jumps (application of operator a on the state vector) 
on the other hand mix the different evolution branches 
and thus tend to reduce the entanglement. 

The fullness of this very intricate stochastic dynamics 
can be captured only by a simulation performed on the 
product basis of the mode-particle Hilbert space, which 
leads to high dimensionality even for a single particle 
[ill . IT^ . Fortunately to understand the central physics, 
such a brute force approach is not necessary and we 
will present a systematic method for identifying a much 
lower dimensional but sufficient subspace in this immense 
Hilbert space. It is spanned by separable particle-field 
states, which represent the full quantum trajectories very 
well after an initial transitional period. The idea should 
be generally applicable in situations where two quantum 
mechanical systems interact and one of them is dissipa- 
tive — a situation ubiquitous in cavity QED, but also 
e.g. in a multi-level atom where the dissipative internal 
degree of freedom and the motion arc coupled by a spa- 
tially dependent pump. 

As the cavity field is well represented by a stand- 
ing wave cosine mode with wave number f{Kx) — 
cos^{Kx), the Hamiltonian commutes with the parity op- 
erator so that the dynamics does not mix the symmet- 
ric and antisymmetric subspaces of the particle Hilbert 
space and we will restrict ourselves to the symmetric part 
only. To further simplify the mathematics in the follow- 
ing we add a sufficiently strong classical part to the co- 
sine potential, so that with Vq,Uo < we are allowed 
a harmonic approximation f{Kx) « 1 — (Kx)'^. At the 
same time we redefine the detuning Ac — Uq Ac and 
V{X)=\Vo\ + \Uo\ X. 

Let us now construct a subspace spanned by separa- 
ble states, which works well in the regime of moderate 
coupling in the sense that the expectation value of the 
projector to the subspace is close to unity for most of the 
time on a trajectory. As a particle with a fixed state vec- 
tor cannot get entangled with the field, we get a coherent 
steady state in the mode. Therefore, if we fix the parti- 
cle in a harmonic oscillator eigenstate, the corresponding 
stationary field is a coherent state |a), which can be sim- 
ply calculated as the particle merely causes a frequency 
shift of the mode by \Uo\ K"^ {x"^) . This coherent state 
in the mode creates an additional average potential for 

the particle reading V ^|a|^^ (Kx)"^. For the combined 

potential we then can again calculate the corresponding 
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FIG. 2: (Color online) System time evolution using an en- 
semble of 300 trajectories started from a product state of vac- 
uum for the mode and (a) |0,^o), (b) (|0,^o) + |2,Co)) /V^. 
After a transitional period during which the photon num- 
ber grows from to its quasi steady-state value the sys- 
tem arrives into the subspace £^'^°''', where it stays from 
then on, so that ^pt^^ii)^ is close to one. Parameters: 

(«;,Vo) = (10,-100)t^,ec, (Ac,7?,f/o) = (0, 2.5, -10)k. This 
Uo value represents a moderate coupling regime for our pur- 
poses. -PoT*"' represent the projectors to only the states corre- 
sponding to m = 0, 2. As we see, the population in these two 
states in the quasi steady-state differs from that in the initial 
condition (50 — 50%): the m — state has higher population, 
in which the cavity cooling effect fis!] is manifested at work. 



particle eigenstate. Iterating this process (or simply solv- 
ing a corresponding nonlinear equation for a) we find the 
subspace 



£(™'^)=span{|a„)|m,e|„,„|2)} 



mGN 

where the particle state is a harmonic oscillator eigen- 
state with oscillator length ^\a^\2 determined together 
with am self consistently. 

In Fig. [5] we demonstrate that gC'^oh) jg ^ very suitable 
subspace to represent the full coupled particle-field dy- 
namics in the moderate coupling regime. This is done 
by calculating the expectation value of the subspace's 
projector P for a state obtained from Monte Carlo wave- 
function simulations of the full coupled dynamics with a 
dimension of « 2000 

Being based on self-consistent product states, the 
above method for obtaining contains essential 

parts of the nonlinear field dynamics but completely dis- 
regards the expected particle- field entanglement. Hence 
this ansatz (and the subspace £(™'^)) is unsuitable for 
stronger coupling where we have to use an extended 
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approach involving a somewhat larger set of separable 
states. While the particle states will still be harmonic- 
oscillator eigenstates in self-consistent average potentials, 
the corresponding mode states have to be determined by 
a new Master equation based on an effective nonlinear 
Hamiltonian invoking all the powers of the photon num- 
ber operator. Note that although it is eliminated from 
the dynamics it is the very quantum nature of the mo- 
tion of the linearly polarizable particle, which creates this 
nonlinearity. Besides giving a very good intuition to the 
underlying dynamics, it is a significant boon of this ap- 
proach that we are left with a dynamics to be solved 
solely on the mode Hilbert space, with a much smaller 
dimension. We start the derivation by casting the Hamil- 
tonian ^ into the form 

Hho = ^ (a^a) ^fet^ + -Ac a^a + irj (a^ - a) . (2) 

The first two terms of that Hamiltonian in the harmonic- 
oscillator approximation can be conveniently expressed 
using the ladder operator defined as & = (^"to^^+^^ata® 
p)/\/2, so that it is an operator on the complete mode- 
particle Hilbert space. Note for future reference that b, 
still obey the usual commutation relations. Recall that 
the frequency ^l{a''a) is an operator on the mode Hilbert 
space. 

Consider now the following basis: \n, m) = \n) |m, 
where \n) is the nth Fock state of the mode, while |m,^„) 
is again the mth harmonic oscillator eigenstate corre- 
sponding to the oscillator length determined by the state 
of the mode (c.f. also Fig. [1]). Due to this dependence 
the basis is not a direct product of bases in the subsys- 
tems' Hilbert spaces. Nevertheless, it is an orthogonal 
basis, since the Fock states for different n are orthog- 
onal, while for a given n the |to,^„) states are orthog- 
onal for different m since they have the same oscilla- 
tor length. We now partition the Hilbert space to sub- 
spaces with a given m £,„ = {|n-,™)}„gN with projectors 
Pm = J2n l'^'"^) {n,m\- 

Let us project the Master equation 37l(i/Ho) to 
the subspace £„■ It is an invariant subspace of the 
first and second terms of the Hamiltonian 

(in fact, the \n,m) states are eigenstates), so that 
e.g. = X;„F,„i?('"')F™. This is not true, how- 

ever, for the pump term (third term of the Hamilto- 
nian ^) and the Liouvillean, since operator a mixes 
the subspaces with different m. Indeed, a |n, m) oc 
|n-l)|m,^„), and (to', ^„_i | to, ^ for aU m', 
therefore {n — l,TO'|a|n, m) 7^ 0. In a first approxima- 
tion we neglect the coupling between the subspaces 
by the operator a and assume H w PmHPm for the 
complete Hamiltonian ^ and the Liouvillean. A mathe- 
matical justification of this is that ^„ is only slowly vary- 
ing with n, so that the corresponding eigenfunctions are 
almost orthogonal. All approximations are a posteriori 
confirmed by numerical tests. Within this approximation 
the dynamics can be solved separately in the different 
Em subspaces. Note that the projected field operator is 
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FIG. 3: (Color online) Wigner functions calculated from the 
stationary field density matrix n'^^^^ . for diflterent vibrational 
particle states m. Same parameters as for Fig. [2] but Ac is 
always adjusted to resonance, i.e. in the Hamiltonian (|3]) 
the term linear in a) a vanish. While for m = the state is 
almost coherent the Wigner function becomes more and more 
banana-shaped with increasing m. 



approximated PmaPm J2m V" + 1 \n, to) {n -I- I, to|, 
and that this approximation is not crucial, but it yields 
that the usual a operator will appear in the projected 
Master equation SUl {Hm)- Within this framework we can 
find the steady-state density operator in the subspace 
£„i, which has the dimension of only the field Hilbert 

space, of the form p^^'^ = E„,n' Pm-i,™' ™) ("-'' "^|, 

where p^m-'^ n' formally a density matrix only on the 
mode Hilbert space. It is the steady-state solution of 
(Hm) with effective Hamiltonian 

Hm = ^uj,,ci\Vo\ + \Uo\a^a) (2m + I) 

—Ac a^a + irj (a^ — a) . (3) 

Note that for Wrcc = we recover the standard single 
mode Hamiltonian linear in a^a, which yields pure co- 
herent states in steady state. 

From a Hamiltonian ([T]) linear in a^a, but containing 
the particle-mode interaction, we have hence arrived at 
the conditional Hamiltonian ([3]) nonlinear in a^a operat- 
ing only on the field Hilbert space. In Fig. [3] we illus- 
trate the nonclassical states this Hamiltonian generates. 
As striking example one obtains banana-shaped Wigner 
functions for higher to indices. 

To find a subspace which represents p^'^ well, we adopt 
an idea from the density matrix renormalization group 
method [l4|. We diagonalize the density operator and 
take a (usually small) number of eigenvectors correspond- 
ing to the leading eigenvalues — in fact, we introduce 
a cutoff at a given value e ^ I in the spectrum. Let 
us remark here that for parameters we studied in prac- 
tice, p^'^ was always very close to a pure state (largest 
eigenvalue very close to unity). We obtain the states 
I'^m.i) = J2n ^m.i-.n \n,m), whcrc i indexes the different 
eigenvectors. It runs between I and Nm, which latter 
is a number depending on e and it may be different for 
different to indices. 

Our most salient approximation is that at this point 
we approximate these states by separable states: 

|*m,i) ~ |0m,i) |"^,C(0„,,|ata|0,„,,>) = |em,i) , 

with \(j>m,i) = ^ ^7n,i-n |«) ■ (4) 
n 
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then slowly leaks out. The particle is started from the 
ground state wave packet of the complete potential, so 
that the system is initially in a product state. Note that 
via the generated entanglement the particle is capable to 
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FIG. 4: (Color online) Projector expectation values in the 
strong coupling regime (Ac, i], Uo) = (0, 2.5, — 100)«:. Ensem- 
ble of 300 trajectories. With e = 10~^ = 1 for all m, 
while with e = 10"" iVo = 3, iV2 = 2, Nm>2 = 1. In both 
cases (P) > ^P'"^"*^'^, but indeed in the second case (P) « 1 

after a transitional period. In thin lines results from a typical 
individual trajectory are displayed to show that the quality 
of £'™''' fluctuates wildly on a single trajectory, while that of 
£ remains close to unity even in this case. 



Note that the particle state has been taken to be an mth 
harmonic oscillator eigenstate corresponding to the oscil- 
lator length calculated from the average photon number 
of the mode state \4'm,i)- This approximation is on one 
hand again justified by the aforementioned mathematical 
argument, but also by a physical argument: during the 
dynamics the pump term and the quantum jumps mix 
the levels of different n, and a particle with finite mass 
cannot follow this mixing immediately. In addition it is 
a posteriori justified by simulations. For the sake of fur- 
ther intriguing the reader we note that the states \em,i) 
are not even in the subspace to which we projected 
our Master equation in the first place. 
We readily arrive at the subspace 



£ _ span{|e,„,i)}^gj^ .^^^^^^^^ 

spanned by separable states, for which the mode states 
are leading eigenvectors of a density matrix, which is a 
steady-state solution of 9}l(Hm) containing the square 
root of V (at a) and hence all the powers of the photon 
number operator. In Fig.|4]we demonstrate that this is a 
very high quality subspace even in the (for our purposes 
strong coupling) regime where the formerly derived sub- 
space breaks down. As expected, the quality of £ 
can be increased by decreasing the cutoff parameter e. 

In Fig. Owe demonstrate one more experimentally very 
easily accessible situation, which clearly exhibits the non- 
classicality of the field generated by particle-field entan- 
glement. Here to even more distill the central physical 
effect, the cavity is chosen to be unpumped, but a few- 
photon coherent state is prepared in the mode, which 




FIG. 5: (Color online) Time evolution (300 trajectories) from 
an initial product state. Parameters: (k, Vb) = (0.1, — 60)tJrec, 
(Ac,7?,?7o) = (0,0, — 100)k; — the particle is chosen to be 
very light in comparison to previous cases, so that it can be 
influenced more easily by the field. The arising field squeezing 
is measured by — log (As), where As is the smaller eigenvalue 
of the field quadratures' correlation matrix. The particle-field 
entanglement is measured by the negativity of the density 
operator's partial transpose [111. [T5|. 



transform the initial coherent field state into a squeezed 
state. 

We demonstrated that few photon nonlinear optics can 
be implemented even by help of a linearly polarizable 
medium, if one includes spatial dynamics of the medium. 
In the low temperature limit where the medium has gen- 
uine quantum properties, a coherent state input field then 
can get entangled with the motional states. Surprisingly 
this complex coupled dynamics can be described by an 
effective nonlinear field Hamiltonian which can be tai- 
lored to experimental needs. As it operates only on a 
restricted Hilbert space it can be efficiently simulated. 
In our example we could identify an 0(10) dimensional 
subspace of the complete 0(1000) dimensional Hilbert 
space to which the system dynamics was confined with 
high probability. While we have concentrated on a sin- 
gle quantum particle here, an analogous dynamics will 
emerge in the many particle case if one invokes only the 
lowest few collective excitations of the ensemble. This 
case will also show a strongly increased nonlinearity via 
collective enhancement and open the route to nonlinear 
optics with single photons. 
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